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The  macroscopic  quantum  theory  of  the  electromagnetic  field  in  a  dielectric  medium  interacting 
with  a  dense  collection  of  embedded  two-level  atoms  fails  to  reproduce  a  result  that  is  obtained 
from  an  application  of  the  classical  Lorentz  local-field  condition.  Specifically,  macroscopic  quantum 
electrodynamics  predicts  that  the  Lorentz  redshift  of  the  resonance  frequency  of  the  atoms  will  be 
enhanced  by  a  factor  of  the  refractive  index  n  of  the  host  medium.  However,  an  enhancement  factor 
of  (n2  -T  2) / 3  is  derived  using  the  Bloembergen  procedure  in  which  the  classical  Lorentz  local- field 
condition  is  applied  to  the  optical  Bloch  equations.  Both  derivations  are  short  and  uncomplicated 
and  are  based  on  well-established  physical  theories,  yet  lead  to  contradictory  results.  Microscopic 
quantum  electrodynamics  confirms  the  classical  local-field-based  results.  Then  the  application  of 
macroscopic  quantum  electrodynamic  theory  to  embedded  atoms  is  proved  false  by  a  specific  example 
in  which  both  the  correspondence  principle  and  microscopic  theory  of  quantum  electrodynamics  are 
violated. 
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I.  INTRODUCTION 


The  effect  of  a  dielectric  host  medium  on  the  spon¬ 
taneous  emission  rate  of  a  two-level  atom  remains  an 
interesting  and  challenging  problem  in  quantum  optics 
with  an  importance  that  befits  a  rigorous  test  of  our  un¬ 
derstanding  of  the  interaction  of  light  with  matter.  The 
essential  characteristic  of  the  problem  is  its  multi-scale 
nature  with  the  atom  being  a  creature  of  microscopic 
quantum  electrodynamics  and  the  dielectric  manifesting 
in  the  realm  of  classical  continuum  electrodynamics.  In 
principle,  it  is  possible  to  represent  both  the  atom  and 
dielectric  microscopically  although,  in  practice,  the  spon¬ 
taneous  emission  rate  is  calculated  using  a  macroscopic 
quantum  electrodynamic  theory  EMSili  in  which 
the  continuous  dielectric  medium  is  incorporated  into 
a  medium-assisted  electromagnetic  field.  While  the  di¬ 
electric  renormalization  of  the  spontaneous  emission  rate 
of  an  embedded  atom  generates  considerable  interest,  a 
mere  handful  of  papers  discuss  the  effect  of  a  dielectric  on 
the  Lorentz  redshift  of  the  resonance  frequency  of  two- 
level  atoms.  Knoester  and  Mukamel  [7j  used  a  Hopfield 
[§]  polariton  model  of  macroscopic  quantum  electrody¬ 
namics  and  found  that  the  Lorentz  redshift  is  enhanced 
by  a  factor  of  the  refractive  index  n  compared  to  the 
vacuum  value.  Crenshaw  and  Bowden  Q  derived  the  en¬ 
hancement  factor  (n2  +  2) /3  using  the  Bloembergen  [l(| 
procedure  in  which  the  classical  Lorentz  local-field  con¬ 
dition  is  applied  to  the  optical  Bloch  equations.  The  two 
procedures  also  produce  different  values  for  the  renor¬ 
malization  of  the  field  (Rabi  frequency)  that  drives  the 
dynamics  of  atoms  in  a  dielectric. 

In  this  paper,  we  derive  the  generalized  optical  Bloch 
equations  for  a  dense  collection  of  two-level  atoms  in  a 
dielectric  host  material  using,  first,  the  Bloembergen  pro¬ 
cedure  based  on  the  classical  Lorentz  local-field  condition 
and,  second,  Ginzburg  macroscopic  quantum  electrody¬ 


namics.  The  derivations  are  simple  and  direct  and  based 
on  well-established  physical  theories,  yet  produce  con¬ 
tradictory  results  for  (i)  the  Lorentz  redshift  of  the  reso¬ 
nance  frequency  of  two-level  atoms  densely  embedded  in 
the  dielectric  and  (ii)  the  renormalization  of  the  field  that 
drives  the  dynamics  of  the  embedded  atoms.  The  micro¬ 
scopic  quantum  electrodynamic  procedure  described  by 
Crenshaw  and  Bowden  11]  confirms  the  Bloembergen- 
based  results.  Then,  we  must  conclude  that  the  macro¬ 
scopic  quantum  electrodynamic  theory  applied  to  dielec¬ 
trically  embedded  atoms  violates  both  the  correspon¬ 
dence  principle  and  microscopic  quantum  electrodynam¬ 
ics. 


II.  DIELECTRIC  LOCAL-FIELD  EFFECTS 


In  the  view  of  Lorentz,  classical  continuum  electrody¬ 
namics  is  better  expressed  in  terms  of  an  atomistic  model 
of  discrete  particles  embedded  in  the  vacuum  that  in¬ 
teract  with  the  microscopic  electromagnetic  field  at  the 
point  of  the  particle  [lj,  113]  •  The  local  field  that  acts  on 
that  particle, 


E  l  =  E  + 


(2.1) 


is  comprised  of  the  macroscopic  Maxwell  field  E  and  the 
reaction  field  of  all  other  particles,  expressed  in  terms  of 
the  macroscopic  polarization  P.  For  a  linearly  polariz¬ 
able  material,  the  polarization  P  =  pNy E^  is  the  prod¬ 
uct  of  the  microscopic  polarizability  p,  the  dipole  number 
density  Nd,  and  the  local  field.  Using  the  Lorentz  local- 
field  condition  EU  to  eliminate  the  microscopic  local 
field  produces  the  Clausius-Mossotti-Lorentz-Lorenz  re¬ 
lation 


£  —  1 

7+2 


(2.2) 
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between  the  polarizability  and  the  macroscopic  dielectric 
constant  e  =  1  +  47 tP/E. 

The  local-field  principle  also  applies  to  nonlinear  me¬ 
dia.  Bloembergen  [loj|  investigated  nonlinear  optics  in 
the  presence  of  a  linear  host  medium  and  found  that  a 
local-field  dielectric  enhancement  factor  of 


where  N  is  the  number  density  of  atoms.  Collecting 
terms  in  V  and  using  the  Clausius-Mossotti-Lorentz- 
Lorenz  relation  (12.21)  yields 

V=^£+^2NpR21.  (2.5) 

47T  o 

Substituting  the  polarization  envelope  (12.51)  into  the  gen¬ 
eralized  Bloch  equations  (12.41)  produces 


accompanies  each  appearance  of  a  macroscopic  field  in 
the  nonlinear  susceptibility.  In  our  notation,  e  =  n2 
for  a  dielectric  and  £  refers  to  the  specific  quantity  in 
the  preceding  equation  and  not  to  any  other  constant 
or  variable  representing  a  local-field  factor.  Bowden  and 
co-workers  [T3,  HE  HI  predicted  intrinsic  optical  bista¬ 
bility  in  a  dense  collection  of  vacuum-embedded  two-level 
atoms  due  to  an  inversion-dependent  local-field  shift  of 
the  resonance  frequency.  Later  work  [9]  reported  the  ef¬ 
fect  of  embedding  the  dense  two-level  systems  in  a  host 
dielectric.  The  dynamics  of  the  two-level  systems  are 
described  by  the  optical  Bloch  equations  [13].  Building 
on  Bloembergen’s  work  0,  Bowden  and  Dowling  [l8j 
showed  that  the  field  that  drives  the  atoms  is  the  local 
field  (12.11).  With  that  substitution,  one  obtains  the  gen¬ 
eralized  optical  Bloch  equations  flq. 


dR.21 

9t 


i(u> 


u>o)R 


ip 
~  2 h 


W  - 

(2.4a) 


dW 

~dt 


ip, 

h 


£?21 


-7||  (W-Weq).  (2.4b) 

Here,  fields  are  represented  in  the  plane- wave  limit  by 
envelope  functions  such  that  P  =  |( Ve~lut  +  c.c.), 
E  =  ^(£e~ZUJt  +  c.c.),  and  El  =  \{£L^~lut  +  c.c.).  The 
macroscopic  spatially  averaged  atomic  variables  in  a  ro¬ 
tating  frame  of  reference  are  R21  =  ( P2ielult)sP ,  R12  = 
(pi2e_*“t)sp,  and  W  =  R22  —  R11  =  (p22)sp  —  (pn)sp, 
where  ( •  •  ■  )sp  corresponds  to  a  spatial  average  over  a  vol¬ 
ume  of  the  order  of  a  resonance  wavelength  cubed  and  the 
Pij  are  the  density  matrix  elements  for  a  two-level  sys¬ 
tem  with  a  lower  state  |1)  and  an  upper  state  |2).  Also, 
/i  is  the  matrix  element  of  the  transition  dipole  moment, 
assumed  real,  7_l  is  a  phenomenological  dipole  dephas¬ 
ing  rate,  7||  is  a  phenomenological  population  relaxation 
rate,  and  Weq  is  the  population  difference  at  equilibrium. 

For  a  linearly  polarizable  material,  the  polarization  is 
P  =  pNcL^L-  For  atoms  embedded  in  a  linearly  polar¬ 
izable  material,  the  polarization  is  the  sum  of  the  linear 
and  nonlinear  components.  Substituting  the  local  field 
EH  into  the  linear  component,  we  have  the  polarization 
envelope 


V  =  pNd 


+  27V/xI?2i 


^  =i(u-u 0  -  ^Np2£W^j  R21~1£W-1±R2X, 

(2.6a) 

dW  i 

—  =  [p££*R2i  -  pi£R*2 1]  -  7||(W  -  Weq).  (2.6b) 

For  now,  we  assume  that  £  is  real.  Then,  the  local-field  ef¬ 
fect  of  the  dielectric  is  simply  an  enhancement  of  the  driv¬ 
ing  field  £,  or  Rabi  frequency  p£/h,  and  the  inversion- 
dependent  Lorentz  redshift  by  £  Q.  The  decay  rates 
remain  phenomenological. 

The  optical  Bloch  equations  0  are  the  basic  semiclas- 
sical  equations  of  motion  for  an  isolated  two-level  system 
in  the  vacuum.  The  two-level  system  interacts  with  its 
environment  through  the  local  field.  Using  the  classical 
Lorentz  local-field  condition  shows  the  effects  of  the  envi¬ 
ronment  to  be:  (i)  a  Bowden-Lorentz  redshift  of  the  reso¬ 
nance  frequency  by  (47r /3h)N p2W  due  to  nearby  atoms, 
(ii)  a  Bloembergen  enhancement  of  the  field  £  by  £  due 
to  the  linear  host,  and  (iii)  a  Bloembergen- type  enhance¬ 
ment  of  the  Bowden-Lorentz  redshift  by  £  due  to  the  lin¬ 
ear  host.  Further  consideration  of  the  generalized  Bloch 
equations,  derived  using  the  classical  Lorentz  local-field 
condition,  is  suspended  until  quantum  electrodynamic 
equations  of  motion  have  been  derived.  At  that  time 
Eq.  EH  will  serve  as  the  basis  for  quantum-classical 
correspondence. 


III.  MACROSCOPIC  QUANTUM 
ELECTRODYNAMICS 

Quantum  electrodynamics  can  be  viewed  as  the  quan¬ 
tized  version  of  Lorentzian  electrodynamics  in  which  dis¬ 
crete  quantum  particles  interact  with  the  vacuum  field 
modes.  When  applied  to  dielectrics,  the  practice  has 
been  to  create  a  macroscopic  version  of  quantum  electro¬ 
dynamics  along  the  lines  of  a  quantized  version  of  con¬ 
tinuum  electrodynamics.  The  macroscopic  theory  can 
be  derived  either  by  quantizing  the  classical  Maxwell 
fields  or  by  applying  a  continuum  approximation  to 
the  microscopic  quantum  electrodynamic  Hamiltonian. 
Ginzburg  [3]  pioneered  the  procedure  of  canonical  quan¬ 
tization  of  the  field  in  a  dielectric  and  applied  it  to 
Cherenkov  radiation.  The  macroscopic  quantization  pro¬ 
cedure  was  limited  to  dielectrics  with  negligible  disper¬ 
sion  and  absorption.  Jauch  and  Watson  [2,LL9|  continued 
Ginzburg’s  work  and  extended  macroscopic  quantization 
to  dispersive  dielectrics  0],  while  other  researchers 
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have  treated  absorption  [2l],  [22|  and  nonlinear  dielectrics 
MM-  Knoester  and  Mukamel  Q  and  Huttner,  Baurn- 
berg,  and  Barnett  [[25]  start  with  the  fundamental  micro¬ 
scopic  Hamiltonian  and  transform  from  coordinate  space 
to  wave-number  space  in  the  continuum  approximation 
to  derive  the  macroscopic  Hopfield  Hamiltonian  [13]  ■  Re¬ 
cent  work  [13,  [H,  [H,  [5o|,  [HJ  includes  macroscopic  quan¬ 
tization  of  fields  in  magnetodielectric  media,  including 
left-handed  negative-index  materials. 

The  spontaneous  emission  rate  of  an  atom  in  a  di¬ 
electric  is  readily  derived  by  the  macroscopic  quantum 
electrodynamic  theory.  In  1976,  Nienhuis  and  Alkemade 
[gj  used  a  macroscopic  version  of  Fermi’s  golden  rule  to 
derive  the  dielectrically  enhanced  spontaneous  emission 
rate  with  the  macroscopic  Ginzburg  fields.  Huttner,  Bar¬ 
nett,  and  Loudon  [53].  amongothers,  have  combined  the 
macroscopic  Fano-Hopheld  [§i,  [33}  theory  with  Fermi’s 
golden  rule  to  derive  the  dielectric  renormalization  of  the 
spontaneous  emission  rate  of  an  impurity  atom,  while 
other  studies  begin  with  macroscopic  Green’s  functions 
[51}  or  auxiliary  fields  (35}.  For  the  most  part,  the  col¬ 
lective  attention  is  focused  on  the  dielectric  renormal¬ 
ization  of  the  spontaneous  emission  rate.  Knoester  and 
Mukamel  Q  obtained  the  dielectric  effect  on  both  the 
Lorentz  redshift  and  the  spontaneous  decay  rate  by  de¬ 
riving  operator  equations  of  motion  from  the  Hopfield 
model.  In  this  section,  we  derive  equations  of  motion  us¬ 
ing  Weisskopf-Wigner  theory  applied  to  the  macroscopic 
Hamiltonian  in  terms  of  the  Ginzburg  field  operators  [361] . 
A  term-by-term  comparison  with  the  generalized  Bloch 
equations  that  were  derived  in  the  preceding  section  un¬ 
der  the  Lorentz  local-field  condition  exposes  an  extraor¬ 
dinary  degree  of  disagreement  between  two  known  and 
accepted  treatments  of  the  effect  of  a  dielectric  host  on 
the  electrodynamics  of  two-level  atoms. 

The  principal  product  of  the  macroscopic  quantization 
theory  is  the  medium-assisted  field  operator 


1/n,  which  is  squared,  and  the  D  =  n3  density-of-states 
factor. 

Local-field  effects  of  a  dielectric  are  suppressed  in  the 
macroscopic  quantization  procedure  and  such  local-field 
effects  must  be  introduced  phenomenologically.  The 
paradigm  that  emerged  from  the  propagation  studies  of 
Hopf  and  Scully  [37]  and  Bloembergen’s  work  in  non¬ 
linear  optics  is  that  the  effect  of  a  dielectric  host  is  to 
multiply  each  occurrence  of  the  dipole  moment  of  a  two- 
level  atom  by  a  local-field  enhancement  factor  0,  So¬ 
using  the  Lorentz  virtual  cavity  model  of  the  local  field, 
the  spontaneous  emission  rate 


T  =  n£2 


3  c3h 


=  n 


(3.3) 


for  atoms  in  a  dielectric  scales  as  n5  for  large  n.  For  an 
atom  in  a  real  cavity,  the  local-field  enhancement  factor 
is  based  on  the  Onsager  model,  and  the  modified  spon¬ 
taneous  emission  rate 


F  =  n 


3  n2 

2  n2  +  1 


2 

To 


(3.4) 


scales  as  n.  A  study  of  local-field  effects  by  de  Vries 
and  Lagendijk  ||38]  found  that  the  Lorentz  virtual-cavity 
model  is  appropriate  if  the  atom  goes  into  a  crystal 
substitutionally  but  that  the  Onsager  real-cavity  model 
should  be  used  for  interstitial  impurities. 

The  macroscopic  quantum  electrodynamic  theory  can 
be  used  to  derive  additional  consequences  of  the  dielectric 
host  for  two-level  atoms.  Taking  the  field  in  a  coherent 
state,  the  effective  Hamiltonian  is 


He ff  =  ^2  +  ft  '}  U lajai 


js 


IX 


ih 


££( 

js  IX 


S?V+o,e<k‘-r*  -gfa\ale-^ 


where  a\  and  a  1  are  the  macroscopic  creation  and  destruc¬ 
tion  operators  for  the  field  modes  and  u>i  is  the  frequency 
of  the  field  in  the  mode  l.  Also,  V  is  the  quantization 
volume,  ek,A  is  a  unit  vector  in  the  direction  of  the  po¬ 
larization,  and  A  denotes  the  state  of  polarization.  The 
spontaneous  emission  rate  of  an  impurity  atom  in  a  di¬ 
electric  can  then  be  obtained  by  applying  Fermi’s  golden 
rule, 

0-7 r 

r=T\(f\Hint\i)\2D,  (3.2) 

to  the  effective  interaction  Hamiltonian  Hjnt  =  —  jia  •  E, 
where  i  labels  the  initial  state  and  /  denotes  all  avail¬ 
able  final  states.  As  is  typically  calculated,  the  dielectric 
renormalization  of  the  vacuum  spontaneous  emission  rate 
of  an  atom  is  found  to  be  n  [j,  0,  0]  due  to  the  dielec¬ 
tric  renormalization  of  the  electric  field  operator  (EH)  by 


(3.5) 

For  species  a,  a3  is  the  inversion  operator  and  aJ±  are 
the  raising  and  lowering  operators  for  the  jth  atom,  g /  = 
(2ltUl  / hvy^fla  (x.  j  ■  e^x)  is  the  coupling  between  the 
atom  at  position  rj  and  the  radiation  field,  xy  is  a  unit 
vector  in  the  direction  of  the  dipole  moment  at  rj,  oja  is 
the  transition  frequency,  and  ga  is  the  matrix  element  of 
the  transition  dipole  moment. 

Except  for  coefficients  of  1/n,  and  macroscopic  field¬ 
mode  operators,  the  effective  Hamiltonian  is  the  same  as 
the  microscopic  Hamiltonian  for  identical  two-level  atoms 
in  the  vacuum.  Equations  of  motion  can  then  be  de¬ 
rived  in  the  same  manner,  which  is  the  primary  reason  for 
adopting  the  macroscopic  formalism.  The  formal  integral 
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of  the  Heisenberg  equation  of  motion  for  the  field-nrode 
operators  is  used  to  eliminate  these  operators  from  the 
remaining  Heisenberg  equations  of  motion.  One  obtains 

=  -iUaoi(t)  +  Timt)e 

1  IX 

n  do  .  ,  . 

IX  j ,s 

+  4E^  sCO  f  dt'e-^-OgfaKt')  (3.6) 
n  ix 

and  a  similar  equation  of  motion  for  the  inversion  op¬ 
erator.  The  procedure  to  evaluate  these  terms  for  two- 
level  atoms  in  the  vacuum  is  generally  known  and  will 
be  considered  in  detail  in  the  next  section.  For  now  it  is 
sufficient  to  note  that  the  only  differences  from  the  vac¬ 
uum  case  are  the  coefficients  of  powers  of  n  and  the  n3 
renormalization  of  the  density-of-states  in  a  dielectric  in 
which  the  sum  over  modes  is  evaluated  as 

f0  dwiu?  J dn{}  (3.7) 

in  the  mode  continuum  limit.  Transforming  to  a  rotating 
frame  of  reference  and  performing  a  local  spatial  average, 
one  obtains  the  Bloch-like  equations  of  motion 

cH?2i  (  47t  ,  \  iu„  - 

—  =i(ur-u.-  -N^„W)  R2i  -  W£IV 

-nX^-R2 1,  (3.8a) 

r)W  2  —  — 

—  =  [iia£*R2i  -  lialRl  1]  -  nT0(W  +  1),  (3.8b) 

where  W  =  ( a3)sp  and  i?2i  =  {-ia-)sp. 

Equations  (13.81)  are  the  equations  of  motion  for  a  two- 
level  atom  in  a  dielectric  host  medium  that  are  derived 
using  macroscopic  quantum  electrodynamics.  The  effects 
of  the  host  appear  as  an  enhancement  of  the  decay  rates 
and  the  Bowden-Lorentz  redshift  by  a  factor  of  the  re¬ 
fractive  index  n,  when  compared  to  the  vacuum  n  =  1 
case.  However,  the  equations  of  motion  derived  using 
the  classical  Lorentz  local-field  condition,  Eqs.  (EH,  dis¬ 
play  an  enhancement  of  both  the  Bowden-Lorentz  red- 
shift  and  the  field  by  a  factor  of  t. 

The  fact  that  local-field  effects  are  suppressed  in  the 
macroscopic  theory  is  well  known.  The  accepted  practice 
10,113  is  to  phenomenologically  associate  a  local-field 


factor  of  £  with  each  occurrence  of  the  dipole  moment  /1 
based  on  the  propagation  studies  of  Hopf  and  Scully  [37| 
and  Bloembergen’s  work  11a  in  nonlinear  optics.  In  this 
case,  Eqs.  (13.81).  become 

=  ,  (a*  -  -  %NAnew)  ft,  -  ^l£W 

—  n£2^I?2i,  (3.9a) 

^  =  -r  [hJ£*R2  1  -  gat£R*21]  -  nPv0(W  +  1). 

(3.9b) 

The  ad  hoc  local-field  correction  gives  the  Bloembergen 
enhancement  of  the  field  that  was  derived  classically  in 
Eq.  (EH,  although  the  redshift  is  over-corrected.  More 
significantly,  the  Bowden-Lorentz  redshift  retains  the  ex¬ 
traneous  factor  of  the  refractive  index  due  to  the  macro- 
scopically  quantized  fields.  Then  the  equations  of  motion 
for  two-level  atoms  in  a  dielectric  host,  Eqs.  (I3.8[>  and 
(EH,  that  were  derived  using  the  macroscopically  quan¬ 
tized  fields  are  inconsistent  with  the  generalized  Bloch 
equations  (12.61)  of  the  preceding  section,  and  the  appli¬ 
cation  of  quantized  macroscopic  fields  to  the  electrody¬ 
namics  of  two-level  atoms  is  contraindicated  by  classical 
Lorentz  local-field  theory. 

IV.  MICROSCOPIC  QUANTUM 
ELECTRODYNAMICS 

In  the  preceding  two  sections,  we  derived  generalized 
Bloch  equations  of  motion  for  two-level  atoms  in  a  di¬ 
electric  host  by  two  well-known  methods  and  obtained 
contradictory  results.  The  most  fundamental  theoretical 
approach  is  to  represent  both  the  atoms  and  the  dielec¬ 
tric  microscopically  0,  0,  m,  E3,  0,  EO,  EH,  0, 0, 0 

and  derive  equations  of  motion  from  first  principles.  We 
show  that  the  results  of  microscopic  quantum  electro¬ 
dynamics  affirm  the  Lorentz  local-field  theory  with  re¬ 
spect  to  the  field  renormalization  and  the  Lorentz  red¬ 
shift.  The  macroscopic  quantum  electrodynamic  theory 
produces  different  results  for  these  effects  and  is  therefore 
not  valid. 

A  dielectric  host  containing  one  or  more  two-level 
atoms  is  modeled  quantum  electromagnetically  as  a  mix¬ 
ture  of  two  species  of  atoms,  a  and  b ,  embedded  in  the 
vacuum.  To  emphasize  the  symmetry  of  the  local-field 
interaction,  both  species  of  atoms  are  initially  treated  as 
two-level  systems.  Species  b  is  later  taken  in  the  har¬ 
monic  oscillator  limit  that  is  associated  with  a  relatively 
large  detuning  from  resonance.  Then,  the  total  Hamilto¬ 
nian  is  comprised  of  the  Hamiltonians  for  the  free  atoms, 
the  free-space  quantized  radiation  field,  and  the  inter¬ 
action  of  the  two-level  systems  with  the  free-space  quan¬ 
tized  electromagnetic  field.  The  multipolar  and  minimal- 
coupling  Hamiltonians  are  related  by  a  canonical  trans¬ 
formation  and  either  can  be  used.  However,  due  to 
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the  canonical  transformation,  the  circumstances  of  the 
rotating-wave  approximation  (RWA)  are  different  for  the 
two  Hamiltonians  [H] .  In  the  typical  derivation  of  the  di¬ 
electric  susceptibility  from  the  minimal-coupling  Hamil¬ 
tonian,  the  RWA  is  invoked  implicitly  by  replacing  po- 
lariton  eigenenergies  with  photon  energies  0,  IH] .  We 
take  the  direct  route  and  use  the  multipolar  Hamilto¬ 
nian  in  the  RWA.  Using  a  plane-wave  expansion  of  the 
electromagnetic  field 

E  =  H  c ■)  Sk.A.  (4-1) 

ix  v 

the  multipolar  RWA  Hamiltonian  is  [HE! 

js  ns  IX 

js  IX 


-  (frr«+°le<k,'rB  -  hraU-e-^)  ,  (4.2) 

ns  IX 


where  aj  and  ai  are  the  creation  and  destruction  opera¬ 
tors  for  the  field  modes  and  uii  is  the  frequency  of  the  field 
in  the  mode  l  with  ki  =  u>i/ c.  For  species  a,  <j33  is  the  in¬ 
version  operator  and  are  the  raising  and  lowering  op¬ 
erators  for  the  jth  atom,  g3  =  (2i tlui/HV)1/2  fia(x-j  •  e^y) 
is  the  coupling  between  the  atom  at  position  r j  and  the 
radiation  field,  Xy  is  a  unit  vector  in  the  direction  of 
the  dipole  moment  at  r j,  oja  is  the  transition  frequency, 
and  Ha  is  the  dipole  moment.  For  species  b ,  g3,  /i[‘ , 

rn,  u>b,  and  fib,  perform  the  same  functions.  Also,  V  is 
the  quantization  volume,  e^y  is  the  polarization  vector, 
and  A  denotes  the  state  of  polarization.  The  polariza¬ 
tion  indices  on  the  variables  have  been  suppressed  for 
clarity.  In  the  two-level  approximation,  the  transitions 
s  =  Am  S  (—1,  0,  +1)  are  treated  separately  and  the  op¬ 
erators  need  not  carry  a  specific  value  for  the  magnetic 
sublevel  pbj. 

Equations  of  motion  for  the  material  and  field-mode 
operators  are  developed  in  a  straightforward  manner 
from  the  Hamiltonian  using  the  Heisenberg  equation.  We 
have 


dai 

dt 


-tunai  +  Yt  siVe-41"'  +  Y 
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(4.3a) 
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dt 


-iujaa3_  °3 a*elkrrj  > 

ix 


(4.3b) 


d<x  3 
dt 


-2  Y(gl<aiel*rT]  +  gfalole-^)  , 
ix 

(4.3c) 


dt 


-iubS™ 


Yh?$a*eik' 

ix 


(4.3d) 


dcn  /  \ 

-^  =  -2^  +  /ir<4^e-ik‘-r"J  . 

.  (4-3e) 

Bloch-like  operator  equations  of  motion  are  obtained  by 
substituting  the  formal  integral  of  the  field-mode  opera¬ 
tor  equation  (14.3al). 


ai(t)=ai( 0)e~iuit+  f  dt'e 
Jo 


x  Y9fa-(t')e-lkrTj  +5Z/l"*?-(t')e_ikl'r 

\  js  ns 

(4.4) 

into  the  material  operator  equations  of  motion  (14.3bl). 
(I4.3cl).  (I4.3dl).  and  (I4.3ell  [47j.  We  transform  operator 
variables  to  different  rotating  frames  of  reference  in  which 
d3_  =  a3_elUat  and  <f"  =  elUbt  are  slowly  varying  quan¬ 
tities.  Performing  the  indicated  substitution  into  Eq. 
(14.3bD  produces 

r1rrj 

~^  =  Y  9i  °3 (0 eikl'Ti at (0) e_i^‘ ~Ua)t  +  a{(t)  £  g3 
IX  IX 
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dt' 


£sf*i(tVMr'-r‘) 


q_e -i(ub-Ua)t 
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dtf 


xYh?*Z-(t'Wkl<rj~rn)  (4-5) 

ns 

in  normal  ordering  with  r^j  =  r j  —  r j.  Likewise,  one 
obtains 
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xY9i*a-^e 
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ik,-(rn-rm) 


(4.6) 


from  Eq.  (14.3dl).  Equations  of  motion  for  the  inversion 
operators  are  obtained  from  Eqs.  (14.3cl)  and  (|4.3cl)  in  a 
similar  fashion. 

The  field  that  drives  an  atom,  Eq.  (El,  consists  of 
the  vacuum  held,  the  self-held,  and  the  reaction  held  and 
we  can  identify  the  terms  on  the  right-hand  side  of  Eqs. 
(14.51)  and  (14.61)  with  fluctuations  due  to  the  vacuum  held, 
spontaneous  decay  from  the  self-held,  near-dipole-dipole 
interactions  between  same-species  atoms  associated  with 
the  reaction  held,  and  near-dipole-dipole  interactions  of 
an  atom  with  the  atoms  of  the  other  species,  also  asso¬ 
ciated  with  the  reaction  held.  The  usual  procedure  is  to 
limit  consideration  to  only  the  spontaneous  decay  rate 
of  a  single  impurity  atom  by  dropping  the  fluctuations 
and  the  single-species  interactions  for  both  the  bath  and 
impurity  atoms.  These  terms  are  retained  here  because 
they  contain  signihcant  information  about  local-held  ef¬ 
fects  in  dielectrics.  For  concreteness,  we  take  species  b 
to  be  the  bath  atoms  and  species  a  to  be  the  two-level 
impurity  atoms. 


V.  NEAR-DIPOLE-DIPOLE  INTERACTION 


The  near-dipole-dipole  interaction  is  the  basic  mech¬ 
anism  of  the  action  of  the  local  held.  The  Weisskopf- 
Wigner-based  procedure  to  evaluate  the  dipole-dipole 
interaction  for  a  dense  collection  of  identical  two-level 
atoms  was  developed  by  Ben-Aryeh,  Bowden  and  En- 
glund  0,  with  corrections  by  Benedict,  Malyshev,  Tri¬ 
fonov,  and  Zaitsev  [H],  to  investigate  single-species  in¬ 
trinsic  optical  bistability.  The  results  apply  to  both 
species  of  two-level  atoms,  individually,  but  we  work  with 
species  b  in  order  to  maintain  consistent  notation  when 
we  take  the  harmonic  oscillator  limit  of  a  two-level  atom 
and  derive  the  Lorentz  local-held  correction  in  a  dielec¬ 
tric. 

We  consider  a  dense  collection  of  identical  two-level 
atoms  of  species  b  in  which  the  atoms  are  evenly  dis¬ 
tributed  in  the  vacuum  with  a  number  density  Nb .  The 
same-species  interaction 


h 


Then 

IT1{  =  f  d^("»-"0(*-O/ij**£»(t/). 

(A  ‘ 

(5-2) 

Applying  the  typical  Weisskopf-Wigner  procedure  [49|, 
1501.  l51|  in  the  mode  continuum  limit,  one  obtains 


rself 

J1 


2^b\^b\\nu, 

3hc3 


(5.3) 


where 


4^3H2 

3  he3 


(5.4) 


is  the  spontaneous  decay  rate.  For  an  atom  of  species 
6,  initially  in  the  excited  state,  7 b  is  the  spontaneous 
emission  rate  To  into  the  vacuum. 

The  pairwise  interaction  of  atoms  is  carried  in  the  re¬ 
maining  m  ^  n  part  of  the  summation.  In  the  Milonni- 
Knight  [46|  model  of  the  interaction  of  two  identical  two- 
level  atoms,  the  strength  of  the  interaction  depends  on 
the  separation  distance  and  the  magnetic  sublevel  tran¬ 
sition.  Then  |46), 


F\  (R)  =  eiR  (-^ 


i 

B3 


for  A  to  =  ±1  transitions  and 


(5.5) 


w  =  (-f  -  Jr)  <“> 


for  Am  =  0  transitions,  where 


Pb 


3  he3  ’ 


(5.7) 


77  —  kbTnrn  —  bJbTnm  f  C,  kb  —  fo>  | ,  Ynrn  —  Yn  and 

rnm  =  |  rnm  |-  Performing  the  summation  over  the  mag¬ 
netic  sublevels,  the  pairwise  dipole-dipole  interaction  can 
be  written  as 


if  =  ${t)\Pb  Y.  -  rnm/c),  (5.8) 

m^n 


where 


B 


nm 


[(^m  ‘  xn)  (xm  •  nnm)(xn  •  nnm  )}Fpkb 

T  nm) 


xYhr*Z-(t')eikl<rn~rm)  (5-1) 

ms 

can  be  extracted  from  Eq.  (14761) . 

The  self-interaction  of  the  nth  atom  with  its  own  re¬ 
action  field  is  characterized  by  the  term  rm  =  r„  in  the 
interaction  as  a  consequence  of  the  relation  ?3  (t)sH(t)  = 
between  Pauli  spin  operators  for  the  same  atom. 


4“  (^m  *  hnm)(  Xn  *  (kbT  nm  )  (5-9) 

and  n„m  =  r„m/r™,  is  a  unit  vector  in  the  direction  of 
rnm  =  r„  —  rm  [1J,  L48(].  Further,  Bnm  incorporates  a 
view  factor  to  account  for  the  arrangement  of  the  dipoles 
in  the  volume. 

The  atoms  are  evenly  distributed  with  a  number  den¬ 
sity  Nb.  For  the  nth  atom,  the  single-species  dipole- 
dipole  interaction  is  obtained  in  a  summation  over  all 
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other  atoms  of  species  b.  In  the  region  near  rn,  the  inter¬ 
action  is  evaluated  by  taking  the  location  of  dipoles  as 
discrete,  while  the  continuum  approximation  is  applied 
elsewhere.  Then 

m:rnrn<o 


+  ^3  {t)\pbNb  (  BZ-{t  -  |r|/c)d3r,  (5.10) 

z  -Iv -vs 


where  5  is  the  radius  of  a  small  spherical  volume  Vs ,  larger 
than  a  cubic  wavelength,  about  the  point  rn.  For  cu¬ 
bic  symmetry,  the  field  generated  by  the  localized  atoms 


The  atom  n  is  located  at  the  origin  of  a  cylindrical 
volume  of  thickness  L  and  radius  Rq.  The  near-dipole- 
dipole  interaction  is  obtained  by  evaluating  the  integral 


rdd 


-M  ?3"(i)  /  # 


rL/ 2 


dz 


-L/2 


pdpK- 


(5.11) 


excluding  a  volume  (4/3)7 t53  about  the  origin  from  the 
range  of  integration,  resulting  in 

^  m  -47»  ^|2e7fc^?(,)<,_(,-|r|/c))sp.  (5.12) 

In  the  limit  <5  — >■  0,  the  near-dipole-dipole  interaction 

j/d  =  (i)?_  (5.13) 

remains  finite.  Here 

n  =  ^NM2  (5.14) 

is  the  strength  of  the  near-dipole-dipole  interaction  and 
<T_  =  (c_(t  —  |rj/c))sp  represents  a  spatially  averaged 
quantity.  The  details  of  this  calculation  can  be  found 
in  the  articles  by  Ben-Aryeh,  Bowden,  and  Englund  14] 
and  by  Benedict,  Malyshev,  Trifonov  and  Zaitsev  [4a] . 

The  atoms  of  species  b  can  be  treated  as  harmonic 
oscillators  if  all  excitation  frequencies  are  far  from  res¬ 
onance  with  uib-  In  this  limit  the  atom  essentially  re¬ 
mains  in  the  ground  state  such  that  sg  —>  —  1.  Then 
the  near  dipole-dipole  interaction  reduces  to  the  Lorentz 
local-field  correction,  shifting  the  resonance  frequency  by 
47rA/'b|/Xb|2/(3/z).  The  microscopic  result  is  in  full  agree¬ 
ment  with  the  classical  Lorentz  local-field  correction  and 
has  been  experimentally  validated  [E3,  m  by  selective 
reflection  of  Rb  from  a  sapphire  window. 


Finally,  all  of  the  results  of  this  section  can  be  applied 
to  the  other  species  of  atom.  Repeating  for  species  a 
yields 
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3  he3 
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7/d  «  (f  —  |r|/c))sp  =  -ivaaJ3(t)a-.  (5.16) 


In  addition, 


7  a  =  2  f3a 


3  He3 


and 


are  defined  for  later  use. 


(5.17) 


(5.18) 


VI.  INTERSPECIES  INTERACTION 

The  effect  of  the  interspecies  near-dipole-dipole  in¬ 
teraction  can  also  be  evaluated  microscopically  using 
Weisskopf-Wigner  theory.  A  single  rotating  frame  of  ref¬ 
erence  is  used  for  both  species  of  atoms  by  making  the 
transformation  Applying  the  results 

of  the  preceding  section,  we  have 

ds  .,  .  ^  pb  r—  .  .  _  7 b  ~n 

—  =  -*(u»  -  wa)<-  -  -fb  +  ^Vb,-  y?_ 


ix  is 

(6.1) 

in  the  harmonic  oscillator  limit  ^  (t)  — >  —1,  where  the 
fluctuating  field 

fb  =  5^(27 ruih/V)1^  ■  ekyk‘  Tnai( o)e-‘("«-»)t 
ix 


is  associated  with  the  spontaneous  decay  rate  by  the 
Kramers-Kronig  relations. 

The  interspecies  interaction  in  Eq.  EH)  describes  how 
a  specific  host  atom  n  interacts  pairwise  with  each  of  the 
impurity  atoms.  The  formal  integral  of  the  equation  of 
motion  of  the  host  atoms,  Eq.  EH),  is 


e(t)  =  -/  dt'e-^-^^f-  -  f  dt'e 
Jo  ">  Jo 
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where  a  =  uit,  —  u>a  —  Ub  —  i'yb/Z-  Substituting  Eq.  (16.21) 
into  Eq.  m ,  one  obtains 
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atoms  in  the  continuum  limit.  For  nonidentical  atoms, 
the  pairwise  interaction  goes  as  [54,  55] 

l2dA  =  ^  BnitUt" -rni/c) 

£ 

i:rni<o 

+  t:\/ PaPbNa  f  Bd-{t"  —  jr|/c)d3r.  (6.6) 

1  Jv-v5 

Then,  Eq.  (16.61)  is  evaluated  as  in  Sec.  V  to  obtain  0 
IddA  =  -^z-Nan*aHb(a-{t"  -  |r|/c))sp 


x  e-i(Wl,-uia)(t*-t'")  Y^g\*dl{t'"y^"Tni.  (6.4) 

is' 

The  term  containing  the  fluctuations  has  been  dropped 
from  consideration  because  the  procedures  presented  here 
apply  only  to  slowly  varying  quantities  and  because  there 
will  be  no  contribution  from  the  random  fluctuations  af¬ 
ter  averaging. 

The  direct  dipole-dipole  interactions  between  impu¬ 
rity  atoms  was  derived  in  Section  V.  Equation  (16.41)  con¬ 
tains  two  such  dipole-dipole  interactions  between  non¬ 
identical  atoms  that  are  integrated  over  the  different 
subspaces  corresponding  to  (i)  bath  atoms  and  (ii)  two- 
level  impurity  atoms.  Due  to  the  relation  a^{t)al{t)  = 
— aJ_(t )  between  Pauli  spin  operators,  the  term  i  =  j  is 
the  special  case  that  is  associated  with  the  renormaliza¬ 
tion  of  the  spontaneous  decay  rate.  This  separates  the 
interaction  I2  =  l2dd  +  i|elf  into  I |elf  for  the  case  i  =  j 
and  l2dd  for  the  summation  over  the  rest  of  the  impurity 
atoms.  The  two  parts  of  the  interaction  will  be  consid¬ 
ered  separately. 


A.  Dielectric  Mediated  Dipole— Dipole  Interaction 


The  nth  atom  of  the  dielectric  interacts  pairwise  with 
every  impurity  atom.  The  term 


rddA 
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l'  A' 


x  J2  (6.5) 
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extracted  from  Eq.  (US),  can  be  evaluated  in  the  same 
manner  as  in  Section  V,  except  that  the  atoms  are  of  dif¬ 
ferent  species.  The  summation  represents  the  effect  of  all 
the  impurity  atoms  on  a  single  atom  of  the  host  material. 
The  sum  over  the  i  ^  j  impurity  atoms  is  performed  (i) 
in  the  near  region  by  the  discrete  summation  over  the  im¬ 
purity  atoms  and  (ii)  elsewhere  by  treating  the  impurity 


=  -^Nanluta-.  (6.7) 

The  quantity  ct_  =  ( d-(t "  —  |r|/c))sp  is  slowly  varying  in 
time  and  the  temporal  integral 

IddB(?)  =  £  dfe-^'-^  -)  Nanlnbd-  (6.8) 

can  be  performed  in  the  adiabatic-following  approxima¬ 
tion.  Repeatedly  integrating  Eq.  (16.81)  by  parts  [56j,  the 
series  can  be  truncated  at  the  first  term  in  the  expansion 
if  the  time  rate  of  change  of  d-  is  much  smaller  than  ad- 
yielding 

=  (v)  (tf)  JWf-  <“> 


The  remaining  part  of  the  I$d  integration  is  another 
interspecies  dipole-dipole  interaction.  In  this  case,  the 
summation  imparts  the  effect  of  all  the  atoms  of  the  host 
dielectric,  modified  by  interspecies  interaction  with  the 
impurity  atoms,  on  the  jth  impurity  atom.  Combining 
terms, 

Id  d(t')  =  -^^Nbfilfia-^-Nanluba^d-  (6.10) 

on  a  on 

becomes 

Id  d{t')  =  ,  (6.11) 

where  \b  is  the  linear  susceptibility  of  species  b.  Adding 
the  direct  near-dipole-dipole  interaction  from  Eq.  (ESI), 
we  obtain 


Idd  -  ivaa3{t)a-{t) 


vacr3(t)a-{t) 


n2  +  2 

=  -i — - — vaa3(t)a-{t).  (6.12) 

Comparison  of  Eq.  (16.121)  with  the  single  species  dipole- 
dipole  interaction,  Eq.  (15.161).  shows  that  the  effect  of  the 
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dielectric  host  is  to  enhance  the  interaction  by  a  factor 
of  £  =  ( n 2  +  2)/3.  Taking  the  local  spatial  average,  W  = 
(o’3)sp  and  R21  =  (—ia-)Sp,  one  finds  that  the  Lorentz 
redshift 


n2  +  2  4n 
— 3  3 h 


7V/x2 


(6.13) 


integration  over  the  subspace  of  two-level  atoms.  Like¬ 
wise,  the  summation  over  the  magnetic  sublevels  invokes 
Sss>-  Then, 

=  -^^PaPbCT-  E  |[(xj  ■  n„j)(x„  ■  n nj)]2F%(karnj) 


is  consistent  with  the  Lorentz  local-field  calculation,  Eq. 
(12. Gap,  while  the  macroscopic  quantum  electrodynamic 
result,  Eq.  (13.9al).  is  not. 


B.  Dielectric-Enhanced  Spontaneous  Decay  Rate 

Most  of  the  elements  of  the  microscopic  theory  of  the 
spontaneous  decay  rate  of  an  atom  in  a  dielectric  are 
common  to  the  treatment  of  the  dipole-dipole  interac¬ 
tion.  In  order  to  show  this  clearly,  we  consider  an  equiv¬ 
alent  derivation  of  the  dielectric  mediated  dipole-dipole 
interaction.  Performing  the  temporal  integrations  first, 
the  interspecies  interaction  (16.41)  can  be  written  as 

/2  =  -  e  9i<4v)*8(u>t  -  u>a)hr^i^~Tn) 

nslX 

X(^)  E  ^(^-Wa)E^*5-(t/)e<k"'(rB"ri)- 

'  is'l'X'  is' 

.  (6-14) 

Applying  the  Milonni-Knight  interaction  with  a  view  fac¬ 
tor  results  in 

h  =  — E  { [(*i  '  “nj)(xn  •  n nj)]F2(karnj) 


+  [(Xj  •  x„)  -  (Xj  ■  n„j )  (x„  •  hnj )]  Fi(fcarnj)} 


+  [(xj  •  x„)  -  (xj  •  n„j)(x„  •  hnj))2  F2(karnj)y  (6.16) 

The  microscopic  treatments  of  the  spontaneous  decay 
rate  in  a  dielectric  0  m  are  missing  elements  of  Eq. 
(16.161)  and  can  neither  affirm  nor  contradict  the  macro¬ 
scopic  theory  of  quantum  electrodynamics.  Because  the 
dielectric  renormalization  of  the  spontaneous  decay  rate 
does  not  have  a  classical  local-field  condition-based  ana¬ 
log,  we  do  not  consider  it  further. 


VII.  DIELECTRIC-ENHANCED  FIELD 


The  dielectric  has  an  effect  on  an  applied  electromag¬ 
netic  field  that  can  also  be  evaluated  microscopically. 
Taking  the  field  in  a  coherent  state,  the  partial  Hamilto¬ 
nian  is 


j 


n: 


:<rV(a 


jt  kr) 


ih 

"2 


e-i(upt- kp-r„) 


n 


S7*  <;n  gi("pt-kPT',)j 


(7-1) 

where  uip  is  the  nominal  frequency  of  the  field  and  f la  = 
[ia£/fo  and  fib  =  /ib£/A  are  Rabi  frequencies.  The  total 
Hamiltonian  is  now  comprised  of  the  Hamiltonians  (ED 
and  (14.21).  Developing  Heisenberg  equations  of  motion 
and  eliminating  the  field-mode  operators  and  the  dielec¬ 
tric  operators  results  in  the  appearance  of 


{[(* 


i)(£n 


i)]F2(karni) 


h  =  TjTO- i(t)£(t)  -  E  j  dt 1 

nslX  ® 


[(Xi  •  X„)  -  (Xj  •  n„,;)(x„  •  n„j)]  F1{karni)}a'l_.  (6.15) 


Xhn*ei{^-K 


f  dt"ea{t'-t")^LS(t") 
/  0  2/i 


(7.2) 


Converting  the  sums,  excluding  i  =  j,  to  integrals  and 
integrating  over  the  subspace  of  two-level  atoms  and  then 
over  the  subspace  of  oscillators  is  equivalent  to  the  deriva¬ 
tion  of  the  dielectric  mediated  dipole-dipole  interaction 
that  was  presented  in  the  preceding  subsection. 

The  renormalization  of  the  spontaneous  decay  rate  of 
a  dielectric-embedded  two-level  atom  is  derived  from  the 
interspecies  interaction  m  in  the  same  fashion  by  tak¬ 
ing  the  target  atom  to  be  the  same  as  the  source  atom. 
The  summation  over  the  two-level  atoms  is  evaluated 
with  the  use  of  the  delta- function  Sjj ,  rather  than  the 


as  an  addition  to  Eq.  (16.31).  Equation  ED  contains  the 
same  type  of  interaction  that  was  evaluated  in  the  previ¬ 
ous  section.  Performing  the  adiabatic-following  approxi¬ 
mation  and  the  sum  over  polarizations,  bath  atoms,  and 
magnetic  sublevels  in  the  mode  continuum  limit,  we  ob¬ 
tain 


h 


n2  +  2  na  ■ 


(7.3) 


The  electromagnetic  field  is  enhanced  by  the  same  factor 
of  £  as  the  reaction  field. 
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VIII.  OPTICAL  BLOCH  EQUATIONS  FOR 
EMBEDDED  ATOMS 


The  macroscopic  optical  Bloch  equations  can  be  de¬ 
rived  from  the  quantum  electrodynanric  equations  of  mo¬ 
tion  in  the  limit  of  large  numbers.  Combining  Eqs.  & 
(10. IIP,  and  (17. 2 J).  one  obtains 

=  ~iivaala-  +  £^Se~l{ujp~Ua)ta33,  (8.1) 

neglecting  the  Gaussian  noise  source  with  zero  mean  and 
absorption.  The  equation  of  motion  for  the  inversion 
operator 


dal 

dt 


nvab\d-  -  ^4^e“i(Wp_a,a)t + Hc- 


(8.2) 


is  derived  in  a  similar  manner.  Optical  Bloch  equations  of 
motion  are  obtained  by  transforming  to  a  frame  rotating 
at  the  frequency  of  the  field  and  taking  a  local-spatial 
average,  as  in  Sec.  III.  We  compare  the  optical  Bloch 
equations 

=  *21  ~  ^ 

(8.3a) 

dW  i 

—  =  --[gam?21-^aO;i],  (8.3b) 

that  were  derived  from  first  principles,  to  the  Lorentz 
local-field-based  equations  (USD-  Based  on  a  favorable 
comparison  of  the  local-field  enhancement  of  the  Lorentz 
redshift  and  the  Rabi  frequency  with  the  classically  de¬ 
rived  result,  we  can  reasonably  assert  that  the  micro¬ 
scopic  theory,  unlike  the  macroscopic  quantum  electro¬ 
dynamic  theory,  satisfies  the  correspondence  principle. 

The  microscopic  theory  allows  us  to  consider  the  more 
general  case  of  of  a  complex  local-field  enhancement  fac¬ 
tor.  Separating  the  real  and  imaginary  parts  of  £,  the 
optical  Bloch  equations  can  be  written  as 

dR2i  (  4ir  ,  \ 

—Qt ~  =  i  (  ujp  -  ua  -  —NffjrW)  i?21 


lihe£W-fnN^WR^ 


(8.4a) 


^  =  idat£*R21  -  iial£Rlx\  -4^a|i?2i|2  (8.4b) 

with  i  =  lr  +  iti.  The  microscopic  theory  justifies  the 
use  of  a  complex  refractive  index  in  the  classical  Lorentz 
local-field  condition.  Then  Eqs.  (ED,  with  phenomeno¬ 
logical  damping,  can  be  derived  by  substituting  the  po¬ 
larization  (12.51)  with  complex  n  into  the  generalized  Bloch 
equations  (12. 4P.  The  imaginary  part  of  the  Lorentz  red- 
shift,  derived  in  this  manner,  was  found  to  be  associated 


with  an  intrinsic  cooperative  decay  for  two-level  atoms 
in  an  absorptive  host  Q .  This  result  is  confirmed  by  the 
microscopic  theory. 

The  optical  Bloch  equations  (18.411  for  dielectric- 
embedded  two-level  atoms  are  derived  from  the  mi¬ 
croscopic  description  of  quantum  electrodynamics  using 
vacuum-based  fields  that  are  known  to  satisfy  the  equal¬ 
time  commutation  relations.  Because  the  field-mode  op¬ 
erators  have  been  eliminated,  the  equal-time  commuta¬ 
tion  relations  cannot  be  discussed  in  the  context  of  the 
optical  Bloch  equations  (18.411  or  Heisenberg  equations 
(18.11)  and  (18.21).  Instead,  the  optical  Bloch  equations,  gen¬ 
eralized  for  a  dielectric  host,  must  demonstrate  conserva¬ 
tion  of  probability.  The  total  population  is  IT2+4|i?2i|2- 
Direct  substitution  from  Eqs.  (18.41)  shows  that  the  tem¬ 
poral  derivative  of  this  quantity  is  nil,  as  required,  in  the 
limit  that  absorption  by  the  atoms  and  the  host  dielectric 
can  be  neglected. 


IX.  SUMMARY 

The  interest  in  the  dielectric  renormalization  of  the 
spontaneous  emission  rate  of  an  atom  embedded  in  a  di¬ 
electric  material  has  obscured  the  inconsistencies  in  the 
macroscopic  theory  of  quantum  electrodynamics.  The 
dielectric  renormalization  of  the  Lorentz  redshift  and  the 
Rabi  frequency,  but  not  the  spontaneous  decay  rate,  can 
be  derived  using  the  classical  Lorentz  local-field  con¬ 
dition  providing  an  independent  check  on  the  validity 
of  macroscopic  quantum  electrodynamics.  The  optical 
Bloch  equations  for  a  dense  collection  of  two-level  atoms 
in  a  dielectric  host  medium  were  derived  using  the  clas¬ 
sical  Lorentz  local-field  condition 

_.  ^  471- 

Ez,  =  e  +  — p 

in  Sec.  II,  while  in  Sec.  Ill,  a  different  set  of  optical  Bloch 
equations  were  derived  using  the  macroscopic  quantum 
electrodynanric  theory.  Both  derivations  are  short  and 
uncomplicated  and  are  based  on  well-established  phys¬ 
ical  theories,  yet  lead  to  contradictory  results  for  the 
Lorentz  redshift  and  the  Rabi  frequency.  If  we  assume 
the  validity  of  the  Lorentz  local-field  condition,  then  the 
macroscopic  procedure  is  proven  to  be  incorrect.  Con¬ 
versely,  the  validity  of  the  macroscopic  quantum  elec¬ 
trodynamic  theory  would  imply  that  the  Lorentz  local- 
field  condition  is  incorrect.  One  deciding  factor  is  that 
the  Lorentz  local-field  correction  has  been  validated  ex¬ 
perimentally  [HU,  [53i| ,  while  the  experimental  record  for 
the  macroscopic  quantum  theory  has  been  inconclusive. 
We  applied  the  more  fundamental  microscopic  theory 
of  quantum  electrodynamics  to  the  same  problem  and 
demonstrated  complete  agreement  with  classical  theory. 
The  differences  in  the  Rabi  frequencies  can  be  reconciled 
with  a  phenomenological  local-field  factor  applied  in  the 
macroscopic  case,  providing  the  virtual-cavity  model  is 
used.  However,  no  such  facile  reconciliation  can  be  pro¬ 
vided  for  the  Lorentz  redshift.  We  conclude  that  both 


the  correspondence  principle  and  microscopic  quantum  turn  electrodynamic  theory, 
electrodynamics  are  violated  by  the  macroscopic  quan- 
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